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Abstract 

We investigate the collective response function and the energy- weighted sums (EWS) 
nik for isovector mode in hot nuclei. The approach is based on the collisional kinetic 
theory and takes into consideration the temperature and the relaxation effects. We 
have evaluated the temperature dependence of the adiabatic, Ei = y^mi/m_i, and 
scaling, = y^m^/mi, energy centroids of the isovector giant dipole resonances 
(IVGDR). The centroid energy E3 is significantly influenced by the Fermi surface dis- 
tortion effects and, in contrast to the isoscalar mode, shows much weaker variation 
with temperature. Taking into account a connection between the isovector sound 
mode and the corresponding surface vibrations we have established the j4-dependence 
of the IVGDR centroid energy which is in a good agreement with experimental data. 
We have shown that the enhancement factor for the "model independent" sum rrii is 
only slightly sensitive to the temperature change. 

Keywords: Fermi system, kinetic theory, response function, energy weighted sum, 
isovector giant dipole resonance, relaxation, temperature 
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1 Introduction 

Many features of nuclei are sensitive to nuclear heating. The nuclear heating influences 
strongly the particle distribution near the Fermi surface and reduces the Fermi- surface dis- 
tortion effects on the nuclear collective dynamics Moreover, the heating of the nucleus 
provides the transition from the rare- to frequent interparticle collision regime. One can 
expect that the zero-sound excitation modes which exist in cold nuclei will be transformed 
to the first-sound ones in hot nuclei. Knowledge of the nuclear collective dynamics in hot 
nuclei allows one to understand a number of interesting phenomena, e.g., the temperature 
dependence of the basic characteristics of isovector giant dipole resonance (IVGDR). The 
existence of the IVGDR in the heated nuclei built on the excited states was established a 
long time ago [2]-[Z]- A systematic experimental and theoretical study of the IVGDR in 
hot nuclei provides a significant information on the isovector collective motion at non-zero 
temperatures [HIE]. 
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A good first orientation in a description of tlie IVGDR in liot nuclei is given by a study 
of the isovector response and the relevant energy weighted sums within the quantum RPA- 
like approaches [10J-[14J or the semiclassical kinetic theory [I3]-[TH]. In the present work, 
following the ideology of the kinetic theory, we consider both temperature and relaxation 
effects on the IVGDR characteristics. Our semiclassical kinetic approach ignores both the 
shell and single particle spin effects. Nevertheless, it seems to be quite instructive for an 
investigation of the averaged properties of the many-body systems. In many cases, it allows 
us to obtain analytical results and represent them in a transparent way. There are also 
some conceptual advantages in the use of the kinetic theory. Kinetic approach involves the 
temperature directly into the equations of motion for the distribution function, e.g., the 
temperature is considered here as a dynamic variable. In contrast, in quantum approaches 
the temperature appears after ensemble smearing of the observable quantities and can not be 
attributed to the equation of motion for the wave function. Moreover the kinetic approach 
can be easily generalized to consider the relaxation (damping) processes by introduce of the 
collision integral [IHl [20]. Note that a similar extension of the RPA due to the involving 
of the coupling with 2p — 2h states provides only the fine structure of the giant multipole 
resonance (GMR) [2TI [22] and an additional smearing procedure for the strength function 
has to be used to derive the corresponding coUisional width of the GMR. Note also that the 
non-collisional fragmentation width of the GMR (spreading of the GMR over non-collective 
Ip — Ih excitations in the RPA) does not drive the system toward a thermal equilibrium but 
rather indicates a redistribution of the particle-hole excitations in the vicinity of the collective 
state [13]. In the kinetic approach, this mechanism is presented due to the Landau damping. 

In what follows we use the kinetic approach to study both the temperature and mass 
number dependencies of the averaged characteristics of the IVGDR, such as centroid energies, 
width, isospin symmetry energy, mass coefficient, energy weighted sum (EWS) and EWS 
enhancement factor. The corresponding analysis within RPA and beyond RPA requires a 
large amount of numerical calculations, see e.g. [23j, which not necessary provide a clear 
understanding of above mentioned macroscopic features of the IVGDR. Both the quantum 
RPA and the kinetic approaches use the effective nucleon-nucleon interaction. We apply the 
effective Landau interaction in nuclear interior and the macroscopic boundary condition in 
surface region. The macroscopic boundary condition includes the phenomenological surface 
tension coefficient and thereby substitutes for an effective interaction within the surface 
layer. Note that the effective interaction is usually not quite well-defined near nuclear 
surface because of strong particle density inhomogeneity in this case and the involving of 
the relevant boundary conditions can be used to improve the description. Our goal is also to 
study the conditions for zero- to first- sound transition in presence of the velocity-dependent 
forces and the effect of the thermal Landau damping [15] on the low-energy tail of the 
strength function under different temperature and relaxation conditions. Similar problems 
were considered in Refs. [IS]- [IE] by neglecting the velocity-dependent part of Landau forces 
and using the simplest boundary condition of Steinwedel- Jensen (SJ) model for the wave 
number k = 7r/2i?o, where Ro is the nuclear radius. However, it is well known [21] that such 
kind of boundary condition do not allow the correct description of the ^-dependence of the 
IVGDR energy for light nuclei. A prove of the modified boundary condition for the isovector 
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eigenmode for a finite Fermi-liquid drop plies an important role in our consideration. 

There are different theoretical approaches to describe the temperature behavior of the 
IVGDR width. Below, we will restrict ourselves by the collisional damping and the thermal 
Landau spreading. The alternative approach is the thermal fluctuation model (TFM) in the 
adiabatic coupling scheme [25] which explains the temperature increase of the IVGDR width 
as an effect of the adiabatic coupling of the IVGDR to thermal shape fluctuations. One can 
expect that the combine of both approaches will be able to provide a satisfactory description 
of both the temperature and the mass number dependencies of the IVGDR width [26j . This 
is an object of our forthcoming investigation. 

Finally, we would like to note that we use the IVGDR as an instrument to study the 
isovector motion in a spherical nuclear Fermi-liquid drop. Our claim is to describe the 
general features of collective excitations, such as the A-dependence of the IVGDR energy 
and the isovector sound mode in nuclear Fermi liquid ignoring many quantum effects. In 
particular, our approach provides the possibility to compare the predictions of the standard 
liquid drop model with the Fermi liquid drop one where the dynamic distortions of Fermi 
surface is taking into account. 

In Section 2 we apply the Landau's kinetic theory to the evaluation of the response 
function in a two-component nuclear Fermi liquid. Both the temperature and the relaxation 
phenomena are taken into account. In Section 3 we derive the energy weighted sums as the 
moments of the strength function. Numerical illustrations, summary and conclusions are 
given in Sections 4 and 5. 

2 Response function within the kinetic theory 

To derive the energy- weighted sums for the isovector excitations, we will consider the density- 
density response of two-component nuclear matter to the following external field 

U,,,it) = \oe-'^'q + \y^'r, (1) 
where Aq is the small amplitude, q is the one-body operator 

A A 

j=l i=i 

and C,j is the isotopic index = p (or +1) for proton and ^ = n (or —1) for neutron). The 
response density-density function x{^) is given by [27J 

where the isovector particle density variation Sp_{r, t) = 6p_ = 6pn — Spp is due to the 
external field Uextit) of Eq. ([T]). 



3 



Below we will apply the kinetic theory in {f,p) phase space to the evaluation of the 
response function of Eq. The particle density variation 6p_ is then given by 

Here, 5^ = 2 is the spin degeneracy factor and 6f_{f,p; t) = 6f_ = 5fn — 5fp is the isovector de- 
viation of the distribution function f{r,p; t) from the equilibrium one at certain temperature 
T 



/eq(Cp) 



1 + exp 



T 



-1 



(4) 



where fi and ep = /2m* are respectively the quasiparticle chemical potential and energy, 
m* is the effective nucleonic mass. 

A small isovector variation of the distribution function 5f_ in Eq. ^ can be evaluated 
using the linearized kinetic equation. In the nuclear volume, where inhomogeneity of the 
particle density is small, the quasiparticle concept of the Landau-Fermi-liquid theory [28] 
can be justified. To evaluate 5f_ we will apply the linearized Landau- Vlasov equation, 
completed by a source term dSt{f) for relaxation processes in the following form [191 EQ] 

+ V- Vr-6f_ - V.-(5f/self + f/ext) " Vp^/eq = 5St[/], (5) 

where v = p/m* is the quasiparticle velocity. We point out that the left hand side of 
kinetic equation can be derived by the Wigner transformation from the linearized time 
dependent Hartree-Fock equation in presence of the external field f/ext [29]. The variation 
of the selfconsistent mean field 6Useif in Eq. ([5]) is then given by the Wigner transformation 
to the corresponding mean field of the RPA. The selfconsistent mean field 6Usei{ is related 
to the Skyrme or Landau effective interaction [30]. Within the Landau- Fermi-liquid theory 
the quantity SUgeu can be derived in terms of the Landau's interaction amplitude Vi^t{'P,p') 
as 



/g 
^^^^t'int(p,p')5/_(r,p';t). (6) 

The interaction amplitude fint(P;P') is parameterized in terms of the Landau constants F/ 



as 



^^int {P,P) = ^ F[Pi{p-p'), p = p/p (7) 
z=o 

and the thermally averaged density of states Nt in Eq. ([6]) is introduced to provide the 
dimensionless constants F/ in Eq. ([71). Namely, 

^T = - I M^,^-^. (8) 



with Nq = gpp m* /2 tt^ h^, where pp is the Fermi momentum. 
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The right-hand side 5St[/] of Eq. ([5]) is the Uehhng-Uhlenbeck type colhsion integral 
hnearized in 6f_. The colhsion integral (5St[/] depends on the transition probability of the 
two-nucleon scattering with initial momenta {pi,P2) and final momenta (p'i,p'2)- low 
temperatures T <^ ep, where ep is the Fermi energy, the momenta (^1,^2) and (1^1,1^2) 
localized near the Fermi surface and the relaxation time approximation can be used, see 

Refs. dSlEolEH], 

00 ^ 

lm,l>l ^' 

where ti is the collisional relaxation time and 6f_^im is a component of the /, m multipolarity 
in j>-space of the isovector variation 5f_. Below we will restrict ourselves to the azimuthally 
symmetric case (longitudinal perturbation) where 5f_ depends only on the angle Opq between 
p and g, i.e. Sf_^im is m-independent, see comment after Eq. (fT^ . The partial relaxation 
time Ti in Eq. ([9]), which corresponds to the Fermi-surface distortion of multipolarity I, is 
derived in this CclSG clS 



1 
n 



I 


■00 / 

dep / 


dnpY*{p) 6St[f] 




POO 

1 dep 


j dnpY*ip) 5f_ 



(10) 



Note that there is no term with / = in the sum ([9]) because of the conservation relation 
for the particle number in collision processes. In contrast to the case of isoscalar mode, the 
inclusion of the / = 1 term in the collision integral of Eq. is due to the nonconservation 
of the isovector current, i.e. due to the pn collisions for the counterstreaming neutron and 
proton flows. The numerical analysis shows [ini H?] that the isovector relaxation time ti 
depends only slightly on the multipolarity / > 2 and we will use below the following form 
for the collision integral 

SSt[f] = --5f_\i=, - -5f_\i>2 . (11) 
n T2 

Here the notations / = 1 and / > 2 mean that the perturbation of Sf_\i=i and 6f_\i>2 in 
the collision integral includes only Fermi surface distortions with a multipolarity I = 1 and 
/ > 2, respectively. 

The collisional relaxation time ti in Eq. f llUp is temperature and frequency dependent. 
The temperature dependence of r/ arises from the smeared out behavior of the equilibrium 
distribution function /eq, see Eq. (j4j), near the Fermi momentum fi9[ [28]. The frequency 
dependence of ti is caused by the memory (non-Markovian) effect in the collision integral. 
It can be shown, see Ch. 8 of Ref. [2H], that the presence of fast collective mode changes 
the energy conservation factor in the collision integral 5St[/] and provides the frequency 
dependence of the collisional relaxation time ti. Following Landau's prescription p8], we 
will assume, see also Refs. [HI [HI [32l [33] . 
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The parameter ai in Eq. (fT2|) depends on the A^A^-scattering cross sections. In the case of 
isotropic energy independent cross sections the result for ai and 02 reads |33l El] 



ai = 3e%/4:TT^ hpcqVp 02 = 5 e|/4 vr^ /Ipoq^^F o'av, (13) 

where ei? is the Fermi energy, vp = pp /m* and peq is the bulk density in the nuclear 
interior. The A^A^-scattering cross sections a^v and (T_ in Eq. (fT^ are given by cXav = 
(o"pp + o"nn + 2crpn)/4 and (J- = crpn/2, where cxpp, dnn and o"pn are the cross sections for 
nucleon pairs with relative kinetic energy close to the Fermi energy. The value of ai is 
significantly different for both vacuum and in-medium reduced cross sections. Using the 
vacuum A^A^ cross sections (Tpp = cr^n = 2.5 2.7fm^ and Upn = a^p = 4.8 5.0 fm^, 
one obtains the following estimate of a2,vac = 2.2 2.3 MeV. The vacuum cross section is 
more appropriate in the surface layer of the nucleus. Due to the Pauli blocking effect one 
can expect that the A^A^ cross sections in nuclear matter should be lower than the one in 
free space. Unfortunately, there is a strong uncertainty in the derivation of the in-medium 
reduced A^A^ cross sections [31]. We will use the the following in- medium estimate of Q;2,buik = 
5.4 MeV, see Refs. [IS1E3]. 

In general the partial relaxation time T2 in Eq. ( ITTl) is larger than ri. It is convenient to 
introduce the relation ai = 02/(1 ~ where 1] is the dimensionless parameter. In the case 
of ?7 — > 1 and ai 00, the relative motion of the proton-neutron fluids is not damped. The 
character of damping of the isovector mode depends on the sign of parameter rj. The zero- 
to-first sound transition is only possible for r/ > [IT]. For rj < 0, the relaxation due to Ti 
leads to the faster equilibration of the out phase proton-neutron motion than the transition 
to the first sound. 

At low temperatures only p near the Fermi surface enter Sf_ and the solution of Eq. ([5]) 
can be found in the form 



t) = e^^'^-^^-'\ (14) 



where dfeq/dep is a sharply peaked function at p = pp and i^uj,qip) depends only on the 
direction of p. Moreover we consider the azimuthally symmetric case where i'u},q{p) depends 
only on the angle 9pg between p and q, and expand Vui^ip) in Legendre polynomials as 

00 

l^ui,q{'P) = ^^Pl (cos 9pq) Z// {p) . (15) 

Using Eqs. (ISl), ([61), dHj) and ([I5]), we obtain 



uj -\ I — q ■ V 

T2 



+ Aog-t/= — z/o(p)^o(cos6'pg) + 77— z/i(p)Pi(cos6'pg). (16) 

T2 T2 
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Substituting expressions ([7]) and (fT5|) into Eq. (fT6|) and performing integration in Eq. (fT6!) 
over p, one can come to the following set of equations for the amplitudes i^iip): 

ui{p) + {21 + 1) ^ ^^'^''^"'(^) - ^o(2^ + l)Q'o(;2) 
= 2(2/ + l)7i^o(p)^ [5/0 - Q/o(^)] - ^(2/ + l)r,iPi{p)Qio{z). (17) 



Here, i/; is the averaged amphtude 



and 



1 /• Pi{x)xPi,ix) 1 u; 
(yz/'(-2) = / ' ^ = cosOpq, 7 = , z = s + ^7, s = — . 

2 J z — X T2qv qv 
-1 

For simplicity we will assume 

FU,^o, fUj^o, f;>2 = o. (19) 

Under the condition (fT9l) . the basic equations (ITTll can be solved with respect to the ampli- 
tude uq. After rather simple calculation one can obtain (see Appendix A) 

_ Xin(^, q) . . . 

^0 - -I . TP/ ~ / V^o, l^uj 

1 + ^0 Xin(w, g) 

where the internal response function Xm{^, q) is given by Eq. (lA 6p . 

The amplitude uq is related to the density-density response function of Eq. ([2]). Substi- 
tuting Eq. (fT^ into Eq. ([S]) and using Eq. (fT^ . one obtains 



Using definition ([2]) and Eqs. (1181) and (I2T1) . we obtain the density-density response function 
x(cu, g) for a given momentum transfer q in the following form 

Equation ([22]) (together with (^KM) 

gives a generalization of analogous result of Refs. p^[T7] 
to the case of the velocity dependent (nonlocal) interaction F[ 0. The poles of the response 
function x{^: q) give the eigenfrequencies of collective excitations u = ujji + iujj and satisfy 
the following dispersion relation 

l + F^Xin(^,g) = 0. (23) 
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2.1 Boundary condition 



For finite nuclei, the dispersion relation has to be augmented by the boundary condition. 
The boundary condition can be taken as a condition for the balance of the forces on the free 
nuclear surface 

n- F\g + n- Fs = 0, (24) 

where n is the unit vector in the normal direction to the nuclear surface S, the internal 
force F is associated with the isovector sound wave and Fs is the isovector surface tension 
force. The internal force F is derived by the momentum flux tensor in the nuclear interior 
and can be evaluating directly using the basic kinetic equation ([5]), see Appendix B. The 
isovector surface force Fs is due to the isovector polarization at the nuclear surface 
Both forces n ■ F\g and n ■ Fs in Eq. ([21]) can be represented in terms of isovector shift of 
the nuclear surface and the boundary condition ( 12^ takes the final form of the following 
secular equation the wave number q, see Appendix B, 



C'sym ~l~ 



X 



cq^sym 



3gro(l + katm^ 



(25) 



where x = qRo, Rq = tqA^^^, the parameter hf derives the Fermi-surface distortion effect, 
see Eq. (IB Sp in Appendix B, and Qsym is the isovector surface tension coefficient, see Eq. 
(IB 9p . In the limit Qsym oo, the boundary condition fl^ leads to the boundary condition 
j'i{x) = of the Steinwedel- Jensen model |i44j- We point out that the secular equation (!25ll 
for q has to be solved consistently with the dispersion relation (J23l) for s. 



3 Energy-weighted sums and transport coefficients 



The presence of the nonlocal interaction in Eq. (1221) leads to some important consequences 
for the properties of the EWS mk{q) for isovector mode. Let us introduce the strength 
function per unit volume 

S{uj,q) =lmxiuJ,q)/Tr. (26) 

The EWS are defined by 

CXJ 

rrikiq) = j d{nw) (fkuf S{uj,q). (27) 



Note that the EWS mk{q) for odd k can also be evaluated by use the dynamic polarizability 
Rex(i^, q)- Considering the two limits — * and uo ^ oc one can obtain [11] 

Rex(a;,g)L^o = 2 [m_i(g) + {hujfm.M + •••] (28) 

and 

Rex(c^,g)L^oo = [^M + (M''^3(g) + •.•] • (29) 
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In the case of cold nucleus T = and no relaxation ti,T2 — >■ oo, applying Eqs. (128!) and (129!) 
to the response function fl22l) and using Eqs. (lA 6P - (1A lOp . we recover well-known results 

m 

m'%) = ^. mf\q) = h^^,q^ mf\q) = h^^f^g\ (30) 
Here, peq is the equihbrium particle density, Csym is the isospin symmetry energy 

Csym = fc.ym,vol = ^€^{1 + F^) ^ 60 McV, (31) 

where 6sym,voi is the volume part of symmetry energy in the nuclear mass formula [H], is 
the Fermi energy and m' is the effective mass for isovector mode 

rn* 

'^'^TTf[/3' rn* = m{l + Fj3) (32) 

and the upper index " (0)" indicates that the corresponding quantity is taken for T = with 
Ti,T2 OO. The renormalized symmetry energy C!,y^ in Eq. (!30l) is given by 

C'sym = Csym + 8 e^/15, (33) 

where the last term on the r.h.s. is due to the Fermi surface distortion effect [36] . 

In contrast to the isoscalar mode, the isovector EWS mi{q) of Eq. flHUj) is model depen- 
dent. As can be seen from Eq. flHUj) . the sum mi{q) includes the enhancement factor (for 
nuclear matter) ^7] 

1 + ^^NM = (m/m*)(l + F[/?,), (34) 

which depends on the nonlocal interaction constant F{ 7^ 0. 

In general the inverse EWS m_i derives the static stiffness coefficient. Evaluating the 
distorted wave function |\l/ad) for a static (adiabatic) constrained field ?7ext = Aog + Aq^* (see 
Eq. ([!]) at — > 0), one can evaluate the corresponding variation of the energy 

6E = (^adl H l^ad) - (^cq| H |^eq) = ^T^^Q', (35) 

4 m \_ { 

where H is the non-perturbed Hamiltonian of the nucleus and 5Q is the change of the mean 
value (g) induced by the distorted wave function. Using Eqs. ( !30|) and ( l35l) . we obtain the 
adiabatic stiffness coefficient (per unit volume) Cq^^a in the following form 

_ d'^5E _ 1 

CQ,ad - - ^ - Csy„,/Peq. (36) 

The cubic sum is related to the stiffness coefficient Cq^^c of the scaling approxima- 
tion [38j. Assuming a scaled form of the perturbed ground state wave function \^sc) = 
e~*^^^ |\l/eq)5 one can evaluate the corresponding change of the energy 



(0) 

777. 

5E = (*3c| H |^3c) - (^eql l^eq) = -^^Q^ ■ (37) 
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Using Eqs. (!30l) . the scaled stiffness coefficient Cq^^c takes the following form 



Cc 



m 



(0) 



2m\ 



(0)2 



^sym/ Peq- 



(38) 



Note that the constrained stiffness coefficient Cg^ad of Eq. ( !36l) is not affected by the Fermi- 
surface distortion since the sum m_i contains the static symmetry energy Cgym- This is not 
case for the stiffness coefficient Cq^sc because the renormalized Cgy^^ enters the sum m^, see 
Eq. fl30p . It can be shown ^\ that both the sum and the stiffness coefficient Cq^sc 
depend on the Fermi-surface distortions of multipolarity / < 2. 

The EWS m_i, mi and allow obtaining the adiabatic, i?ad, and scaled, Egc, average 
energies (centroid energies) of IVGDR |40j 



E: 



ad 



m 



(0) 



m 



(0) 

-1 



(0) 



m 



(0) 



m 



(39) 



There is a significant difference between the adiabatic energy, i?ad, i-e., derived by a static 



stiffness coefficient Cgym, and the scaled one. 



Esc, associated with the isovector sound in 



nuclear Fermi liquid. The Fermi-surface distortion effects increase the stiffness coefficient 



and leads to an increase of the centroid energy E^c of the isovector mode. 



The low frequency (see Eq. ( I28l) ) sum m_3 is related to the transport coefficient. Assum- 
ing a slow time dependence of the external field Uextit) and evaluating the corresponding 
time dependent wave function |^I/(t)), one can find the change of the energy in the following 
form 



5E = {^{t)\H\^{t)) 



m 



(0) 



Am 



(0)2 



6Q'. 



(40) 



Using Eq. (HOi) we obtain the transport (mass) coefficient Bq as 



Br 



d^6E 
d 6Q^ 



m 



(0) 



2m 



(0)2- 



(41) 



Using the mass coefficient Bq and the stiffness coefficient Cg^ad (see Eq. f lHBjl ). the eigenfre- 
quency cUmacr and corresponding eigenenergy -Emacr for the macroscopic eigenvibrations can 
be derived as 



Q,ad 



Br 



E„ 



\ 



m 



(0) 



m 



(0)- 



(42) 



We point out that the sum m^g can not be evaluated within the Landau- Vlasov kinetic 



approach by the use of the low frequency expansion of Eq. fl28l) . That is because the Fermi 
liquid gets into Landau-damping regime at u; ~ and special attention must be paid to the 
low energy region in Eq. (|27l) . see below. 
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4 Results and Discussions 



We will present results of numerical calculations based on the consideration of previous 
sections. In this work we adopt the value of tq = 1.2 fm and the effective nucleon mass 
m* is taken as m* = 0.9m which corresponds to the Landau parameter Fi = — 0.3. For 
the isovector interaction parameter Fq we have used Fq = 1.41 to keep a reasonable value 
of isospin symmetry energy Csym of the order of 60 MeV, see Eq. (I3T1) . The interaction 
parameter Fl will be derived and discussed below. For the relaxation parameters in Eq. 
f fT2|) we use the values of 0^2 = 5.4 MeV and 77 = — 0.1 which correspond to the in-medium 
reduced A^A^ cross sections fWi [T7] . For more clear interpretation of some numerical results 
we will also use the relaxation parameters a2 and rj beyond these well-established values. 



4.1 Eigenenergy and enhancement factor in cold nuclei 

The interaction parameter Fl can be estimated by considering the enhancement factor 
1 + ^AfA/ for the isovector EWS see Eq. fl34|) . Following Ref. pRJ], we derive the 

photoabsorption cross section crabs(^^) in terms of the strength function S{uj,q) Eq. fpUj) as 
follows 

, , 471^ NZ , 
CTabsl^^j = — 5 —ujb[uj,q). (43) 

In the case of the velocity independent A^A^-interaction, the cross section crabs('^) PHl) is 
normalized by the ordinary Thomas- Reiche-Kuhn sum rule [13] (see mi[q) in Eq. flHUl) for 

I^NM = 0) 

00 

(0) f Alf. ^ ^ ^ 27r2ne2ArZ „ ^ 

H.TRK = / "-[nw] O-abs(t^) = 1" for I^NM = (44) 







at T = and ri, r2 oo. 

Taking into account the velocity dependence of the A^A^-interaction with Fi 7^ and 
F[ 7^ 0, we note that both the enhancement factor knm 7^ in mi(g) of Eq. flHUl) and the 
corresponding correction at the last term of the boundary condition fIB 17p affect the sum 
rule (jHj). For k^m 7^ 0, using Eqs. fHHj) . flHUj) and fIB 17p . we obtain the following result 
[39] 

~(o) /■ X /X f q[{A)V N Z , 

m\' = I d[nuj) (Tabs(t^) = , — r(l + '^NM) for knm > 0, (45) 

J mc \qo{A) J A 



at T = and ti,T2 00. Here qo{A) and q[{A) are the lowest roots of Eq. fIB 17p for 



i^NM = and K^M 7^ 0, respectively. The value of interaction parameters F[ can be now 
obtained from a fit of the evaluated enhancement factor ^i^V'^i'^trk experimental 
data. In this work we have adopted F[ = 1.1. Our estimate of the enhancement factor is 
about 10% for light nuclei and increases to 20% for heavy nuclei which is in a good agreement 
with experimental data |46j . 
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Fig. 1: Dependence of the IVGDR energy on the mass number A: the dotted hne 1 is the 
first-sound regime without Fermi surface distortion (adiabatic approximation, see i?ad in the 
Eq. (|39i) ): the dotted hne 2 shows the calculation which takes into account the Fermi surface 
distortion up to the multipolarity / = 2 (scaling approximation, see Esc in the Eq. (l39i) ): the 
solid line was obtained by use of the dispersion relation (123|1 and the secular equation (125|) . 
The dashed line shows the result within S J model with commonly used [Til [TBI [T7] value of 
the wave number q = ti/2Rq. The solid dots are the experimental data from Ref. [IHj. 

In finite nuclei, both the IVGDR eigenenergy -Eivgdr and the EWS rrik are the compli- 
cated functions of the mass number A. In contrast to the classical Steinwedel- Jensen model 
[H] , the value £'ivgdr • ^^^^ is not the constant but increases with A [46j . Within our Fermi 
liquid approach, the A-dependence of the IVGDR eigenenergy and EWS occurs due to the 
boundary condition of Eq. fl25l) on the moving nuclear surface. 

In Fig. [T]we show the dependence of the IVGDR energy (multiplied by A^/^) on the mass 
number A. The calculations have been performed for Qsym = 10.5 MeV and F[ = 1.1. The 
solid line is the eigenenergy obtained from the dispersion equation ( |23l) augmented by the 
boundary condition of Eq. fl25|) . Both dotted lines in Fig. [1] have been obtained from the 
EWS definitions of the centroid energies Eq. (IH^ (curve 1 for i?ad and curve 2 for E^c)- The 
significant upward shift of the scaled energy. Esc, with respect to the constrained one, i^ad, 
is due to the Fermi surface distortions of multipolarity / < 2 presented in Egc- An additional 
upward shift of the exact eigenenergy (solid line) is due to the higher multipolarities / > 2 
of the Fermi surface distortions presented in the dispersion equation ( l23l) . As seen in Fig. 
dl we reproduce quite well the average behavior of the IVGDR energy -Eivgdr- Note that 
the boundary condition j[{qRo) = of the Steinwedel- Jensen model [H] as well as the 
commonly used wave number q = 7r/2i?o [El [iHl [TT] do not describe the y4-dependence of 
the IVGDR energy correctly, see dashed line in Fig. [H 
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4.2 First sound regime 

We consider the first sound regime as the displacement of the spherically-symmetric Fermi 
surface without its deformation in momentum space. In the case of the velocity dependent 
effective A^A^-interaction with F{ 7^ 0, the first sound eigenvibrations differ from the cor- 
responding ones in the classical Steinwedel- Jensen model [44j- For the sake of simplicity 
we will consider the longitudinal eigenvibrations assuming [/ext(^) = and no relaxation 
(rfc = 00) in Eq. 1^. Using Eq. ffT^ and expansion (instead of Eq. ffT^ ) 

J^oj,giP) = ^lm{q, Uj)Yi.raip), (46) 

Im 

we will transform the kinetic equation to the following set of equations 

uj i^im - vf G'l, {lm\ q ■ p \l'ni) vvm' = 0, (47) 

I'm' 

where = 1 + F//(2/ + 1) and 

{lm\ q ■ p \l'm') = j dVtp Yi*^{p) cos{9gp)Yi,m'ip)- 

Using Eq. (HM and assuming i'im\i>2 = 0, we obtain from Eq. p7|) the following closed 
equations for amplitudes Uqq and Uiq: 

s 1^00 - -^G[uio = 0, s uio - -^G'qUqo = (48) 

and the corresponding dispersion relation 

u = -^VF qy^G[~G{. (49) 

Finally, using the definition ( 1311) of the symmetry energy Cgym, we obtain the eigenenergy 
energy -Efirst of the IVGDR in the first sound limit in the following form 



E^rst = hh^\i§^q = hJ^il + «:^m) q. (50) 



The energy -Efirst of Eq. ( ]50i) differs from the one E^j of the Steinwedel- Jensen model 



Es, = hJ^q (51) 
V m 

by the enhancement factor 1 + ktva/- Similar result was reported earlier in Ref. 
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Fig. 2: Strength function S{uj,q) from Eqs. ([221) and (^KM for Fi = - 0.3, = 1.1, 
7] = - 0.1, A = 208. Solid line for T = 2 MeV, aa = 5.4 MeV and dashed line for T = 0.5 
MeV, as = 100 MeV. 



4.3 Strength function 

Performing the numerical calculations of the response function x(^) q) f l22l) . one can evaluate 
the strength function S{LJ,q) fl2Ul) and the EWS mk{q) flTH) for finite temperatures T 7^ 
and in presence of relaxation. The strength function S{u},q) is sensitive to the interaction 
parameters and to the relaxation properties. Because of Fq > 0, the IVGDR strength 
function contains both the sound mode contribution at s > 1 and the Landau damping 
region at s < 1. This is illustrated in Fig. O 

To show the presence of the Landau damping in the IVGDR S{uj, q) in a transparent 
way, we have plotted in Fig. [2] the result for the zero-sound regime ujrT2 ^ 1 at T = 0.5 
MeV, ^2 = 100 MeV (dashed line). The Landau damping appears there as a wide bump 
on the left side of the narrow sound peak. For high temperature (solid line in Fig. [2]), the 
sound peak becomes wider due to the decrease of the relaxation time (coUisional relaxation), 
see Eq. (1121) . and due to the collisionless thermal Landau damping, which increases with 
T, see Ref . [T^ . As can be seen from Fig. [21 overlapping of both the sound peak and the 
Landau damping bump leads to the asymmetry of the IVGDR strength function at high 
temperatures. This feature of the IVGDR strength function is observed experimentally [7j. 
Note that the IVGDR width, which can be derived from Fig. [21 represents a collisional part 
of total width only and this one is significantly smaller the experimental width of Ref. [7j. 
Additional part of the IVGDR width is caused by the fragmentation mechanism and we will 
take into account this fact below in Fig. [HI 

Sensitivity of the strength function S{(jJ^ q) to the interaction parameter F[ is demon- 
strated in Fig. 131). The inclusion of the nonlocal interaction F[ 7^ increases the isovector 
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Tico (MeV) 



Fig. 3: Strength function S{uj, q) from Eqs. ^ and (TATOl) for T = 2 MeV, aa = 5.4 MeV, 
Fi = - 0.3, r/ = - 0.1, A = 208. Solid line for F[ = 1.1, and dashed line for F[ = 0. 

stiffness coefficient and shifts the IVGDR to the higher energy. Moreover, since the interac- 
tion parameter Fl enters the enhancement factor 1 + knm of Eq. the photoabsorption 
cross section crabs(i^) P3|) grows with F[ > 0. 

Presence of the 1/ri term in the colhsion integral in Eq. ffTTj) increases the width of 
the IVGDR resonance and does not much affect its energy centroid. However, this term 
influences strongly the zero- to ffist-sound transition for the isovector mode p!7l [18] . In 
general, a decrease of the collisional relaxation time T2 leads to the fast damping of the Fermi 
surface distortions and thereby to the zero- to ffist-sound transition. This is demonstrated 
in Fig. m 

The solid line in Fig. H] shows the numerical result for the photoabsorption cross section 
aahs{Lj) from Eq. fj^Hj) for long relaxation time regime (zero-sound regime, 0^2 = 5.4 MeV) 
at rj = 1, i.e., ti = 00. To show the zero- to ffist-sound transition, we have plotted in Fig. H] 
(dashed line) the cross section crabs(i^) which is obtained in the ffist sound regime ujjiT2 <C 1 
at 0:2 = 0.1 MeV. This transition happens as a shift of the resonance energy to the energy 
of the ffist sound eigenmode, -Egrst, given by Eq. (130]) 

Efirst ^ 17.5 MeV for A = 208. 

This value of -Efirst significantly exceeds the SJ estimate Esj 14.6 MeV for A = 208 
obtained at the boundary condition j[{qRo) = [B], see also Eqs. ( 150|) and ( 15T|) . 

The behavior of the photoabsorption cross section (Jabs(i^) is essentially different for the 
case where the relaxation time T2 exceeds significantly the relaxation time ti. In this case 
the relaxation of the relative proton-neutron motion is faster than the zero- to ffist-sound 
transition and the ffist-sound peak of the IVGDR disappears. This effect is shown in Fig. [5] 
where the dashed line was obtained at 77 = —1, i.e., T2 = 2ri. 
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Fig. 4: Photoabsorption cross section o"abs(t^) from Eqs. (1^ . (HHj) and (lA lOp for T = 2 
MeV, A = 208, r/ = 1, Fi = 0, F{ = 0. Solid line for a2 = 5.4 MeV (zero-sound regime) and 
dashed line for a2 = 0.1 MeV (first-sound regime). 




-hco (MeV) 



Fig. 5: The same as Fig. HJ but for 77 = —1. 
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Fig. 6: Temperature dependence of the "model independent" EWS mi{q) obtained from 
Eqs. f lTTl) . fpUl) and fl^ normalized to the TRK sum rule mi^TRKil)- The calculations 
performed for A = 208 using a2 = 5.4 MeV, 77 = — 0.1 and Fi = — 0.3. Solid line for 
Fl = 1.1 and dashed line for = 0. 

4.4 Energy-weighted sums and centroid energies 

We have studied the temperature behavior of the "model independent" EWS mi{q) and 
the enhancement factor mi{q) / mi^TRKil) , where, see Eq. (|27j) . 

"^i,trk(Q') = J d{fvjj) hhj S{uj,q) for Fi = F[ = 0. 


For non-zero temperatures and in presence of the relaxation, the energy- weighted sum mi (g) 
has been evaluated using the definition fl27l) and the response function xi^^ q) from Eq. (I22p . 
In Fig. [6], we have plotted the ratio mi{q) / mi^TRK{(l) as a function of temperature T. 

As can be seen from Fig. El the enhancement factor is only slightly sensitive to the 
temperature variation. For the nucleus ^°^P6, we have from Fig. [6] the following estimate 

mi(g)/mi,TRK(g) ~ 1-2. 

In Fig. [71 we compare the centroid energies Ei = mijm^x and Ej, = a/ ms/rrii evalu- 
ated using Eq. (1271) . and the IVGDR eigenenergy Eji obtained from Eqs. (!23l) and (!25l) . 

The significant upward shift of energy E^ curve with respect to Ei is due to the Fermi 
surface distortion effect. The cubic sum which enters the expression for E^ is associated 
with the scaling transformation and the quadrupole distortion of the Fermi surface ^2] . On 
the other hand, the distortion of the Fermi surface causes an additional contribution to the 
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Fig. 7: Temperature dependence of the centroid energies Ei = ^J^m^Jrnl'x (solid line) and 
-E3 = \fnv^frn\ (solid line) obtained from Eq. (1271) . and the IVGDR eigenenergy (dashed 
line) obtained from Eqs. ([23]) and (1231) . The calculations were performed for A = 208 using 
a2 = 5.4 MeV, 77 = 2/3, m*/m = 0.9 and Fj' = 1.1. Dotted line shows the classical energy 
E^acr from Eq. (|42l). 

stiffness coefficient (see also Eqs. ([HSD and (IHHl) ) providing a growth of the centroid energy 
£■3. In contrast to E3, the centroid energy Ei is derived by the inverse sum m_i, where the 
contribution from the Fermi-surface distortion effects is negligible, see also Eq. ( 136|) . Note 
that the presence of the Fermi-surface distortion effects increases also the resonance energy 
Er because the dispersion equation (123|) includes all multipolarities of the Fermi surface 
distortion. A small decrease of the energy Eji with growing temperature T in Fig. [71 is due 
to the fact that the Fermi-surface distortion effects become weaker for higher temperatures. 
In contrast to E3, which is nearly temperature independent, the centroid energy Ei increases 
slightly with T. That is because a growth of resonance width with temperature leads to a 
decrease of the inverse sum m_i in agreement with its definition of Eq. ([27|) . 

In Fig. [71 the dotted line shows the behavior of the eigenenergy -Emacr for the macroscopic 
eigenvibrations given by Eq. (112|) . To provide the convergency of the integral in Eq. (1271) 
for the inverse sum m_3 in Eq. ( H2l) . we have introduced the cut off parameter Ecnt in Eq. 
( 127|) . Namely, we have used 

m_3(g) = j d{hLu) (hw)-^ S{uj,q). (52) 

-Bcut 

The cut off energy E^ut was derived from the requirement that the Landau-damping region 
s < 1 should be removed from the EWS since must be related to a given sound 
eigenmode only. Note that the energy interval [i?cut,cxD] must be large enough to provide 
a reasonable exhaustion (> 90%) of sum rule (1^51) . A good agreement of -Emacr with the 
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IVGDR eigenenergy allows one to conclude that the cut off procedure in Eq. (!52|) can also 
be used for the consistent evaluation of the mass coefficient Bq of Eq. fHTl) . 

Finally, it should be noted that there is the limiting temperature T^^ for the IVGDR 
existence. The limiting temperature can be deduced from the decomposition of the Xmi^^, q) 
of Eq. (lA 61) in powers of 1/s for the first sound limit (high temperature regime) at |s| ^ 1. 
The presence of the finite relaxation time ri caused by the collisions for the counterstreaming 
neutron and proton flows leads to the additional temperature dependence of the eigenfre- 
quency which is specific for the isovector mode only. The eigenfrequency of the isovector 
first sound decreases with temperature and disappears at the limiting temperature Tn^. 
The magnitude of the limiting temperature Tum for the IVGDR depends significantly on the 
relaxation parameter rj. For heavy nuclei, the numerical estimate [IT] provides the value of 
Tn^ ~ 7 MeV. 



4.5 Damping and spreading width of IVGDR 

The EWS mk{q) can be used to analyze the spreading of the strength function. We will use 
the following definition of the spreading width [TTl HTj 

mi(g) mo(g) 

7W = TT TT- (53) 

mo{q) m_i(g) 

In the case of a Lorentzian shape of the photoabsorption cross section 

"-^-^"^= (i^-i^o)^ + CT^ ' ^ = ^' 

one can find from Eqs. ( H3i) . ( 1271) and ( l53l) the relationship between the quantities Ei and 
7(g) and the resonance characteristics Eq and F 

Eo = E,, 7(g) = -F(l + 0(F/Eo)) for F/^o « 1- (55) 

TT 

In the case of small damped collective vibrations (F ^ Eq), the resonance collisional 
width F can also be evaluated from the dispersion equation ( 123|) . The solution of this 
equation 

u = uji + iuj, (56) 

defines the energy of the giant resonance Eji = Ttwr and its collisional width Fcoi = —2hwj. 
In the application of the kinetic approach, in particular, of the collision integral 5St[/] to 
the finite nucleus, the difficulty is the derivation of the NN cross sections cTav and (T_ in Eq. 
( |T3l) which become f^-dependent in nuclear interior. The above mentioned (see comment to 
Eq. fll3p ) vacuum and in- medium values of the NN cross sections and thereby relaxation 
parameter a2 in Eqs. f ll2p and 0131) give the lower and upper theoretical limits. The r- 
dependence of the relaxation parameter 02 can be taken into account phenomenologically 
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by introduce of the effective relaxation parameter as following 



"eff 



drpcq{r)/a2{r) 



(57) 



where a{r) is parametrized by 



Ct2(r) = ttvac H TTT ("bulk " ttvac) , 

Peq(O) 

with Ovac = 2.3, MeV and Obuik = 5.4, MeV, see comments after Eq. (fT3|) . The equilibrium 
particle density Peq(r) derived as 



Peq(r) = po 



exp 



where po = (47rrQ/3)~^ and a = 0.6 fm. Below we will replace relaxation parameter a2 in 
Eq. (fT^ by the effective one cteff- 

Besides the collisional width Fcoi, the experimentally observable width of the IVGDR 
includes the fragmentation width. Within the semiclassical kinetic theory, this mechanism 
of resonance spreading can be considered as an additional dissipation due to the single 
particle scattering on the moving surface of the nucleus (one-body dissipation [IHl IMl 150]). 
Instead of r2, we will use the effective relaxation time Teg which contains the contribution 
from both two-body and one-body dissipations. Namely, see also Refs. [T71 [33], 



1 

T"cfr 



1 

T2 



Twall 



(58) 



We will use the one-body relaxation time r^aii in the following form 



Twall 



2-Ro ^ 

- S ■) 

V 



3vf 



^4 - 

6 Kep 



(59) 



The parameter ^ in Eq. fl59|) depends on the model of the one-body dissipation 
SI]. We consider ^ as a free parameter which is determined from a fit of the total IVGDR 
width r to the experimental data at zero temperature T = 0. 

In Fig. [8] we have plotted the temperature dependence of the width F derived from the 
EWS by use of Eq. ( l55l) (solid line 1) and from the dispersion relation applying Eqs. ( l23l) . 
(!56|) . (!57|) and (!58l) (solid line 2). The numerical calculations were performed for the nucleus 
^"^^Sn where the experimental data are known for a wide range of temperatures T = -r 3 
MeV [SIEIE]. 

We point out that considering the experimental data we assume that the nuclear ex- 
citation energy E* is related to the nuclear themperature T by the Fermi-gas formula 
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Fig. 8: Temperature dependence of the IVGDR width F: the sohd hne 1 was obtained from 
the EWS definition given by Eqs. and the sohd hne 2 and the dashed hne were 
obtained from the dispersion equation ( |23l) . The calculations were performed for the nucleus 
^'^^Sn using a^g from Eq. ( 1571) and i] = —0.1. Sohd line 2 for Fi = — 0.3, F[ = 1.1 and 
dashed line for Fi = = 0. 



T = ^y E*/a where a is level density parameter. This fact leads to some uncertainty because 
the derivation of parameter a is model dependent. However, within the kinetic theory, the 
use of the macrocanonical ensemble, and thereby temperature T, is preferable because the 
ensemble smearing is assumed at the derivation of basic collisional kinetic equation. 

In Fig. [HI the total resonance width F grows with temperature mainly due to the temper- 
ature dependence of the collisional relaxation time T2, see Eq. ( JT2l) . As can be seen from Eq. 
f l59|) . the temperature dependence of the one-body dissipation is too weak to be responsible 
for the fast increase of the IVGDR width. The result in Fig. [H] shows also that the collisional 
width does not give a good description of the observable growing of the IVGDR width with 
T. An additional increase of the IVGDR width with temperature can be achieved taking 
into consideration the coupling of the IVGDR to thermal shape fluctuations |25| which do 
not present in our approach. Note that a decrease of the in-medium NN cross section leads 
to an increase of the temperature dependence of the IVGDR width (TUj also. 

For low temperatures, the total resonance width F derived by the EWS (solid line 1 
in Fig. [8]) is close to the one obtained from the dispersion relation (solid line 2 in Fig. [8]). 
That is due to the fact that the evaluated photoabsorption cross section approaches to the 
Lorentzian shape in this case. For higher temperatures, the assumption T/Eq <^ 1 in Eq. 
f l33|) is not fulfilled and the expressions fl3^ and fl3^ for the derivation of F can not be 
correctly used. 
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5 Summary and Conclusions 



A goal of this work is the derivation of the macroscopic characteristics of isovector modes 
in a two-component Fermi-hquid drop from our knowledge of the nuclear IVGDR. Starting 
from the collisional kinetic equation ([5]), we have derived the strength function and the 
energy weighted sums m^, for the isovector excitations in the heated nuclear matter and the 
finite nuclei. An important ingredient of our consideration is the inclusion of the velocity 
dependent A^A^-interaction for both isovector and isoscalar channels. Our consideration is 
valid for an arbitrary collision parameter ai in Eq. (fT2l) and can be used, particularly, for 
the transition region from the zero sound- to first sound (hydrodynamic) regime in nuclear 
Fermi-liquid. 

We have adopted a Fermi liquid drop model with two essential features: (i) The linearized 
kinetic equation is applied to the nuclear interior, where the relatively small oscillations of 
the particle density take place; (ii) The dynamics in the surface layer of the nucleus is 
described by means of the macroscopic boundary condition which is taken as a condition for 
the balance of the forces on the free nuclear surface. 

Due to the consistent solution of the dispersion relation (^Hl) and secular equation 
our model provides a satisfactory description of the A-dependence of the eigenenergy of 
the IVGDR, see Fig. [H In contrast to earlier consideration of Ref. [21] performed within 
the scaling approximation, we take into account all multipolarity / of the Fermi surface 
distortion. Moreover, the value of Landau's interaction parameter F[ k, 1.1 has been derived 
from a fit of the evaluated EWS enhancement factor to the experimental data. Note that we 
do not use a concept of the effective A^A^-interaction within the surface layer of the nucleus 
applying instead the relevant boundary condition and avoiding thereby the uncertainty in 
the derivation of effective interaction in the region of strong particle density inhomogeneity. 

The present study has shown the following: 

1. The Landau damping occurs in the isovector strength function Siu, q) at low temper- 
atures as a wide bump on the left side of the narrow sound peak (see Fig. H]). For high 
temperature, the overlapping of both the sound peak and the Landau damping bump 
leads to the asymmetry of the IVGDR strength function S{uj,q). A similar asym- 
metry of the IVGDR strength function is also observed experimentally at non-zero 
temperatures [7]. 

2. The isovector EWS shows only minor temperature dependence. In particular, the 
"model independent" EWS mi{q) and the corresponding enhancement factor 
Tni{q) /miTRKio) are practically constant in the interval of temperature T = -^ 5 
MeV, see Fig. El 

3. The "model independent" EWS and the enhancement factor mi{q) /miTRKil) 
are slightly sensitive to the relaxation (damping) processes. Note that the correspond- 
ing problem can not be accurately considered within the standard quantum mechanics 
because of non-Hermite Hamiltonian in this case. 



22 



4. The inclusion of the nonlocal interaction F[ ^ increases the isovector stiffness co- 
efficient and shifts the IVGDR energy to the higher values. The lowest order EWS 
m_i(g), mi{q) and m^{q) derive the adiabatic, E'l, and scaled, i?3, energy centroids. 
The centroid energy Ei is close to the classical result of the Steinwedel- Jensen model 
and differs from the first sound limit by the enhancement factor I + ki^m-, see Eq. ( l39l) . 
The Fermi distortion effects do not contribute into the centroid energy Ei. In contrast, 
the scaling energy E^ is associated with the quadrupole distortion of the Fermi surface 
and exceeds significantly the centroid energy Ei. 

5. The often used classical derivation of the IVGDR eigenenergy -Emacr Eq. fH2|) through 
the isovector mass coefficient Bq has to be revised in the case of nuclear Fermi liquid. 
That is because the inverse sum m_3, which enters the mass coefficient Bq, can not 
be directly evaluated within the kinetic theory due to the Landau-damping region at 

~ 0. To provide the convergency of the sum m_3, we have proposed the cut-off 
procedure introducing an appropriate cut off parameter Ecnt into the definition of m_3 
in Eq. ^Tj\ . Due to this procedure we have achieved a good agreement of i^macr with 
the adiabatic centroid energy Ei. In general, such kind of cut-off procedure can also 
be used for the evaluation of the transport coefficients within the Fermi-liquid theory. 

6. The inclusion into the collision integral of the relaxation of / = 1 component (term 
~ I/ti in Eq. f ill I) ) influences strongly the zero- to first-sound transition for the 
isovector mode. In particular, in the case of T2 > ti, the relaxation of the relative 
proton-neutron motion is too fast and the short-relaxation limit a2 ^ does not 
provide the zero- to first-sound transition (see disappearance of the first-sound peak 
in Fig. E]). A similar phenomenon was earlier discussed in Refs. fTli \TE\ where the 
velocity independent Landau's interaction with F[ = has been used. 

7. Our analysis of the IVGDR width performed within the kinetic theory shows that the 
collisional and one-body damping does not reproduce the fast increase of the IVGDR 
width with temperature. The additional mechanisms of damping, e.g., the coupling 
of the IVGDR to thermal shape fluctuations [25] , have to be involved to improve the 
agreement of the temperature dependence of the IVGDR width with the experimental 
data. 

Finally, we would like to note that the semiclassical kinetic approach, used in this article, 
is highly convenient for a study of the averaged properties of the nuclear dynamics. This 
approach provides an information on the macroscopic characteristics without a detailed 
knowledge of the wave function of the nucleus. An important advantage of the kinetic theory 
is that the temperature and the relaxation effects enter the equation of motion directly. 
Many results can be here presented in a clear and transparent form. Here, the claim is 
to describe the general features of collective excitations, such as the A-dependence of the 
IVGDR energy, in a systematic way ignoring many quantum effects, e.g., the shell structure 
effects. In particular, the kinetic approach allows us to compare the results of standard 
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liquid drop model [H] with the Fermi liquid drop one where the dynamic distortions of 
Fermi surface are taken into consideration. 



Appendix A: Internal response function 

Using Eq. (|T7|l . we will evaluate the averaged amplitude z/q- Taking into account Eq. 
( fT9l) . we will reduce Eq. (ITTl) to the following coupled equations 

^ F' ^ 1 

yQ{p)+FQQQQ{z)vQ+^QiQ{z)vi-\QQQQ{z) = i-fuoip)- [1 - Qooiz)]-iri'yiyi{p)Qoo{z), (A 1) 

o Z 

ui{p) + 3FQQio{z)iyo + F[Qn{z)ui-3\oQio{z) = -3i-fUo{p)Qoo{z) -3ivi^iip)Qwiz), (A 2) 
where we have used Qio{z) = z Qoo{z). Solving Eqs. flA ip and flA 2p we obtain 

/ N ^ Xo(^) ^ c./ z Xojz) F{ zxo{z){z-i7] 7) ^ 

sD(z) + ^7X0(2:) sD{z) + i-fXo{z) 3 + 27X0(2:) 

and 

3z s Xo(2:) 3z s Xoiz) ^ 

sD[z) + 27x0(2;) sD[z) + 27X0(2;) 

_ p/ [z'^Xoi.z) + 1/3] [sD(2:) + ^7X0(2;)] -i-j z {z-iri 7)^0(2;) - . . 

D{z)[sD{z)+i^Xo{z)] ^ ^ 

where Xo(2;) = Qoq{,z), 

D{z) = 1 + 32777^X0(2;) 

and the relation Qii{z) = z^xoi^) + 1/3 has been used. Multiplying both Eqs. ( lA 31) and 
f lA 4p by — [5f/(27r/j.)'^A^T] (dfeq/dep), integrating over p and using Eqs. ([H]) and ffTSj) . we 
obtain two closed equations for amplitudes i?o and Pi. Solving then these equations with 
respect to uq, we obtain 

Xin(^, q) , P-v 

Here Xin(a;,g) is the internal response function 

Xin C^, g = Xm,0^, q) - ^1 1 I F/~(2)r ^' ^ 

1 + Fix^^H^>g) 



where 



V ^ / ^/cq(gp) 2:X0(2;) 

^ '^^ NtJ (27rn)3 Sep .D(^) + 27Xo(^) ' ^ ^ 
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and 



1 
1 



gdp dfe, 



z [z 



7) Xo{z) 



gdp df, 



[z'^Xq{z) + 1/3] [sD{z) + z7Xo(^)] - ^7 z{z - ir] 7)Xo(^) 
L>(z) [sD{z)+i-iXo{z)] 



(A 9) 



(A 10) 



Appendix B: Boundary condition 

In this appendix we are going to determine the boundary condition from the balance of 
the forces on the free nuclear surface given by Eq. fl24p . The internal force F in Eq. ([21]) is 
related to the momentum flux tensor 11^13 in the nuclear interior 



(B 1) 



The momentum flux tensor can be evaluated using the basic kinetic equation Taking 
the 1-st p-moment of Eq. ([5]) one can obtain the following expression for the momentum 
flux tensor lic^p [HI 112] 



n, 



al3 



6P 6 



60, 



a/3! 



where 5P is the pressure caused by the iso vector sound wave 



5P 



i / (1^^^ ^f-ifj] t) + ^p,q (5p_(f, t) = Csym 5p_(r, t) 



and Saai3 is the pressure tensor due to the Fermi surface distortion effect 



6a, 



af3 



Here p = {p^^^n + Pcq,p)/2, 
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Vf q 



and X is the displacement field related to the velocity field u. Namely, 



d_ 

dt 



X{r,t) = - {1 + K)u{r, t) , M(f , t) 



gdp p 



Peq J {l-nhf 



m 



Using Eqs. fIFT]) - flBl]) . we obtain 



n-F\ 



Csym Peq ~ gPi^ ) V • X + 2pF^(^ " X) 



9_ 



J r=Ro 



(B 2) 



(B 3) 



(B 4) 



[B 5) 



(B 6) 



(B 7) 
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To evaluate the isovector surface tension force Fs we will consider the variation 6Es,sym 
of the surface symmetry energy caused by the isovector polarization at the nuclear surface 



ds e 



(B 8) 



where Qsym is the coefficient related to the volume, 6sym,voi, and surface, 6sym,surf, terms 
entering to the symmetry energy -Egym in the mass formula 



Q 



9 

•-' sym,vol 



sym 



E. 



sym 



1 {N-zy 

2 A 



sym, vol "sym, surf 



[B 9) 



8 ^sym,surf 

In Eq. (IB 8\} . the parameter ^ is the dynamic isovector shift of neutron-proton spheres in 
units of ro 



(B 10) 



where 



5R,{t)=Roasmio{r)- (B 11) 

The amplitude as{t) of the isovector shift of the nuclear surface is connected with the 
displacement field x- To establish this connection we note that, for the case of sharp nuclear 
surface, the isovector displacement field x is given by [44J 



X(r,t) = ai{t)-V4ji{qr)Y,o{r)] . 

The boundary condition for the normal component of the velocity field reads 

d 



n ■ u\s 



dt 



6Ri{t). 



Using Eqs. (lB6ll . (IFTT!) and (lBT3l) 



x(l + k) ' 



X = qRq. 



The variation 6Es,sym of the surface energy derives the surface pressure 



SPs 



d 6E. 



Pcq 



d SRi 5S 3 ro 
Taking into account Eqs. flFTTll . flFlil) and flBTSl) 

8 Peqflix) 



n-F. = -6P. 



Qsym 5Ri{t). 



Qsymai{t)Yio{r) 



(B 12) 



(B 13) 



(B 14) 



(B15) 



(B 16) 



3gro(l + k) 

Inserting Eqs. fIB 7\i and fIB 16p into Eq. flM|) and using Eq. fIB 12p , we obtain the following 
secular equation 

1 



^CsymPeq r^r'r ' ^2 



2 2 

Pf H — :^Pf 



Pcq 



2 4 

-PF+-- , X 

X 3gro(l + Kj 



Q 



sym 



3[{x) = 0. (B 17) 
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